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and 

B.  HUTCHISON 
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An  iterative-type  thermal  model  applicable  to  forest  canopies  was  tested  with  data  from  two  diverse  forest  types.  The 
model  framework  consists  of  a  system  of  steady-state  energy  budget  equations  describing  the  interactions  of  short-  and 
long-wave  radiation  within  three  horizontally  infinite  canopy  layers.  A  state-space  formulation  of  the  energy  dynamics 
within  the  canopy  is  used  which  permits  a  factorization  of  canopy  geometrical  parameters  from  canopy  optical  and 
thermal  coefficients  as  well  as  environmental  driving  variables.  Two  sets  of  data  characterizing  a  coniferous 
(Douglas-fir)  and  deciduous  (oak-hickory)  canopy  were  collected  to  evaluate  the  thermal  model.  The  results  show  that 
the  model  approximates  measured  mean  canopy  temperatures  to  within  2°C  for  relatively  clear  weather  conditions  and 
deviates  by  a  maximum  of  3°C  for  very  hazy  or  foggy  conditions. 


Introduction 


Rapid  and  accurate  assessment  of 
renewable  resources  is  an  increasingly  im¬ 
portant  task  facing  remote  sensing  spe¬ 
cialists.  Mathematical  abstraction  of  en¬ 
ergy  processes  of  vegetation  canopies  is  a 
useful  technique  for  relating  sensor  re¬ 
sponse  to  environment-canopy  interac¬ 
tions.  Such  an  understanding  is  required 
in  order  to  make  timely  inferences  about 
the  condition  of  forestry  and  agricultural 
resources  from  remote  sensors.  In  the 
thermal  regime,  several  authors  have  re- 
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ported  success  in  estimating  evapo- 
transpiration  of  crops  from  thermal  sensor 
data  (Heilman  et  al.,  1976;  Reginato 
et  al.,  1976;  and  Soer,  1980).  Several 
models  have  been  reported  that  describe 
the  energy  balance  of  vegetation  either  in 
terms  of  a  single  leaf  (Gates,  1968; 
Wiebelt  and  Henderson,  1977 ;  and  Kimes 
et  al.,  1978),  or  an  abstract  layered  canopy 
(Alderfer  and  Gates,  1971;  and  Deardorff, 
1978).  Few  models  have  been  described 
that  characterize  the  energy  flows  within 
vegetation  canopies  as  a  function  of 
the  canopy  geometry  (Goudriaan,  1977; 
Norman,  1979;  and  Kimes  et  al.,  1981). 

0034-4257/81 /040311  +  16$02.50 
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In  this  study,  a  modification  of  the 
thermal  canopy  model  reported  by  Kimes 
et  al.  (1981)  that  is  applicable  to  forest 
canopies  is  described  and  evaluated.  The 
model  is  an  iterative  type  formulation  of  a 
system  of  steady-state  energy  budget 
equations  describing  a  canopy  as  three 
horizontal  layers.  Each  layer  is  described 
in  terms  of  reflective  and  thermal  radi¬ 
ation  coefficients  and  leaf  inclination  an¬ 
gle  distributions. 

The  model  emphasizes  the  radiative 
energy  processes  occurring  both  within 
and  above  the  vegetation  canopy  layers 
and  relates  these  flux  transfers  to  detailed 
consideration  of  the  canopy  geometry. 
Other  energy  transfer  processes  such  as 
sensible  heat  and  evapotranspiration  are 
included  but  in  a  fairly  simplistic  manner. 
Standard  expressions  from  the  literature 
are  utilized.  The  modular  form  of  the 
model  makes  it  fairly  easy  to  replace  the 
existing  expressions  as  warranted.  Thus, 
the  model  is  capable  of  evaluating  partic¬ 
ularly  the  radiative  energy  processes  from 
widely  diverse  canopies. 

There  were  two  broad  objectives  of  the 
work  reported  here.  First  was  to  de¬ 
termine  if  we  could  re-express  a  previ¬ 
ously  developed  thermal  exitance  model 
(Kimes  et  al.,  1981)  in  a  more  usable 
form  which  would  permit  a  wide  variety 
of  engineering-type  target /background 
studies  to  be  performed  in  a  more  practi¬ 
cal  computational  fashion.  The  second  ob¬ 
jective  was  to  evaluate  the  model  with 
validation  data  collected  from  two  differ¬ 
ent  forest  types;  a  Douglas  fir  (Pseudo- 
tsuga  menziesii)  canopy  and  a  mixed  de¬ 
ciduous  canopy. 

Here,  we  first  describe  the  updated 
model  structure  and  solution  .approach. 
This  is  followed  by  a  description  of  the 
experimental  sites  and  methods  for  ob¬ 
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taining  the  required  model  input  data. 
Finally,  the  results  and  summary  are 
given. 

Model  Structure 

The  updated  model  formulation  given 
below  differs  in  three  major  aspects  from 
the  original  model  developed  by  Kimes  et 
al.  (1981).  First,  the  model  has  been  re¬ 
written  in  a  vector- matrix  form  with  a 
state-space  characterization.  That  is,  there 
is  a  specific  identification  of  a  state  vector 
X,  the  canopy  layer  temperatures,  a  con¬ 
trol  vector  U,  the  meteorological  driving 
variables,  and  a  parameter  vector  P.  In 
this  formulation,  the  long-wave  energy 
budget  terms  have  been  factored  into  a 
geometric-dependent  term,  the  S  matrix, 
and  the  long-wave  thermal  source  terms. 
This  factorization  permits  the  precalcula¬ 
tion  of  the  S  matrix  for  a  wide  variety  of 
canopy  situations  and  then  the  subse¬ 
quent  convolving  of  these  matrices  with 
the  U  and  P  vectors  to  evaluate  canopy 
thermal  variations.  A  second  potential  ad¬ 
vantage  of  this  factorization  is  that  it 
results  in  a  linear  system  of  equations 
with  respect  to  S  suggesting  the  use  of 
linear  filtering  theory  to  estimate  S  from 
X  (Sorenson,  1966). 

The  second  modification  is  the  use  of 
more  simplifying  assumptions  in  the 
model,  particularly  with  regard  to  the 
short-wave  absorption  calculation  as  dis¬ 
cussed  later.  Finally,  in  solving  the  non¬ 
linear  energy  budget  equations,  explicit 
use  was  made  of  the  closed  form  expres¬ 
sions  for  the  Jacobian  of  the  system  in  a 
Newton- Raphson  technique. 

In  the  material  that  follows,  individual 
expressions  for  the  component  energy 
budget  processes  are  summarized  and  an 
explicit  expression  for  the  elements  of  the 
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Jacobian  matrix  are  given.  The  geometri¬ 
cal  factorization  of  the  energy  budget 
equation  for  the  long-wave  flux  transfers 
is  derived. 

Energy  balance  framework 

The  model  is  a  plane-parallel  abstrac¬ 
tion  of  a  vegetation  canopy  divided  into 
three  horizontal  layers.  Two  additional 
source  layers  are  given  by  the  atmosphere 
above  the  canopy  and  by  the  underlying 
ground.  An  energy  balance  framework, 
assuming  steady-state  conditions,  is  for¬ 
mulated  for  each  of  the  three  vegetation 
layers  (sinks)  as  a  function  of  the  five 
source  layers.  In  the  expressions  that  fol¬ 
low,  i— 1,2,3  represents  the  sink  or  vege¬ 
tation  layers  and  /  =  1,2, 3, 4, 5  represents 
respectively  the  atmosphere,  the  three 
vegetation  layers,  and  the  ground  source 
layers  of  energy  flux.  The  combination  of 
the  t,  i  indices  thus  represents  flux  from 
source  layer  /  to  sink  layer  i. 

The  vector  expression  for  the  energy 
balance  equations  for  layers  1,2,3  consid¬ 
ering  long-wave  transfers,  short-wave 
transfers,  sensible  heat,  and  evapo- 
transpiration  can  be  written  as 

F(X,P,U)=0,  (1) 

where 

XA(X,X2X3)r=the  average  layer  tem¬ 
perature  vector  for  layers  1,  2,  3; 

^A(6j»  i  =  1,2,3;  alt  i  =  1,2,3;  eg;  R,;  S; 
A)  =  the  parameter  vector  char¬ 
acterizing  the  canopy  layers; 

e(,  Oj  =emissivity  and  absorptivity  of 
the  vegetation  layers; 

eg,ag  =  emissivity  and  absorptivity  of 
the  ground  layer; 

R,  =leaf  stomatal  resistance  to  water 
vapor  diffusion  (Note  that  the  resis¬ 


tance  of  the  boundary  layer  to  diffu¬ 
sion  is  estimated  from  other  mea¬ 
sured  parameters); 

S  =  long -wave  flux  transfer  matrix 
calculated  from  geometrical  proper¬ 
ties  of  the  canopy; 

A = short-wave  flux  absorption  coeffi- 

U*(TaTgWS  RH  S W)T  is  the  control  or 
input  vector; 

Ta  =air  temperature  (°C); 

Tg  =  ground  temperature  (°C); 
WS=wind  speed  (m/sec); 

RH=  relative  humidity; 
SW=short-wave  flux  (w/m2). 

F  may  be  rewritten  in  the  following  ma¬ 
trix  form  which  explicitly  separates  the 
geometrical  properties  of  the  canopy  S 
from  the  remaining  energy  terms  (Smith 
-et  al.,  1981). 

F  ao  B(X)  rS  —  o  B(X) 

+  A  +  H(X)  +  LE(X)=0, 

(2) 

where 

B  =  vector  of  long-wave  emission  terms, 
H  =  vector  of  sensible  heat, 

LE  =  vector  of  evapotranspiration  terms, 
a  =  Stef  an- Boltzmann  constant. 

The  significance  of  this  factorization  is 
that  a  wide  variety  of  abstract  or  cannoni- 
cal  canopies  may  be  characterized  by  pre¬ 
calculation  of  the  S  matrix.  This  matrix 
table  may  then  be  convolved  with  the 
appropriate  meteorological  driving  varia¬ 
bles  to  simulate  diurnal  behavior  for  a 
wide  spectrum  of  scenarios. 

The  vector  equation  (2)  may  be  ex¬ 
panded  for  each  layer  and  the  explicit 
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dependence  on  parameters  or  input  varia¬ 
bles  may  be  indicated  by 

MBiTjSn+BiXjS* 

+  B(X2)S13+B(X3)S14  +  B(rg)S15} 

+Ax-aB(Xl)+H(XliWS,Tl,) 

+LE(XliWS,Ta,Rl,RH)=0, 

(3a) 

tafiBiTjSu+BiXjS* 

+B(X2)S23+B(X3)S24+B(T4)S25} 
+A2  — oB(X2)  +  H(X2;  WS,Ta) 
+LE(X2-,WS,Ta,R„RH)=  0, 

(3b) 

^a3a{B(ra)S31  +B(X1)S32 

+B(X2)S33+B(X3)S34+B(rg)S35} 
+A3-oB(X3)-hH(X3;WS,ro) 
+LE{X3;WS,Ta,Rl,  RH)=0. 

(3c) 

The  formulation  for  each  energy  budget 
component  used  in  the  model  is  given  by 

Long-wave:  B(Xt)=ei(Xi  +  273)4,  (4) 

Qn  nrt-wa  vp 

absorption:  A,  =  ABS(i)-SW,  (5) 

where 

ABS(i)= short-wave  absorption  coeffi¬ 
cient  calculated  by  an  optical 
absorption  model  which  uses 
a  Monte  Carlo  technique  to 
include  multiple  scattering  ef¬ 
fects  (Kimes  and  Smith,  1980; 
Kimes  et  al.,  1980). 

Sensible  heat: 

H{Xi;WS,Ta)=hc{WS)(Xi  —Ta)  (0) 


where  hc,  the  convection  coefficient,  is 
taken  from  Tibbals  et  al.  (1964) 

Evapotranspiration: 

L£(X(;  WS,ra,  B,  RH) 

=  p^[s(X,)-,WS(T„)]MX,), 

(7) 

where 

h(Xt)= latent  heat  of  vaporization  of 
water  at  temperature  X(; 
s(X,)  =  water  vapor  density  inside  the 
leaf  at  saturation  at  the  leaf  tem¬ 
perature  X,  (g  cm-3); 
s(Ta)  —  water  vapor  density  at  saturation 
of  the  free  air  beyond  the 
boundary  layer  of  the  leaf  at  the 
air  temperature  Ta  (g  cm-3); 

Ra  —  resistance  of  the  boundary  layer 
to  water  vapor  diffusion; 

The  terms  s(X(),  s(Ta),  and  Ra  are  calcu¬ 
lated  from  measured  driving  variables,  U 
(Kimes  et  al.,  1981). 

The  discussion  of  the  S  matrix,  which 
controls  the  interception  of  long-wave  flux 
within  the  canopy  layers,  is  given  later. 

Explicit  evaluation  of  the  Jacobian 

An  iterative  Newton-Raphson  Tech¬ 
nique  (Burden  et  al.,  1978)  was  used  to 
solve  the  system  of  nonlinear  thermal 
equations  (2)  since  the  Jacobian  of  the 
system  may  be  analytically  derived  in 
closed  form.  This  method  involves  iter¬ 
ative  evaluation  of  the  following  expres¬ 
sion  about  an  initial  guess  X0  until  8X 
converges; 

8Xi(X-Xo)=0rJ)',Jr[-F(Xo)]. 

(8) 
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where 

J=system  Jacobian  =  (3F/3X]x=Xo. 
The  Jacobian  of  the  system  is  given  by: 


3X„  ’ 

=  2«nemSnm/(Xw+273)3+fin 
x{4em/(Xm=273)3+hJffl 


3f,(Xn) 


+  ^[S(Xn)--RHS(Ta))-± 
h(Xn)  3 s(Xn) 


+ 


R,  +  R<t  3X, 


n,  m  =  l,2,3. 

(9) 


Snm  is  the  Dirac  delta  function.  Given 
specific  algebraic  expressions  for  s(Xn) 
and  h(Xrl)  the  partial  derivatives  are  easily 
evaluated. 


Geometrical  factorization— S  matrix 

A  significant  simplification  of  the  ther¬ 
mal  model  from  the  development  re¬ 
ported  earlier  (Kimes  et  al.,  1981)  is  the 
factorization  of  the  geometric-dependent 
terms  from  the  energy  related  source 
terms  for  the  long-wave  flux  transfer 
processes.  This  factorization  is  made  pos¬ 
sible  essentially  because  of  the  lack  of 
multiple  scattering  in  the  thermal  regime 
between  canopy  components  whose  emis- 
sivities  (absorptivities)  are  assumed  nearly 
unity  and  by  the  fact  that  the  thermal 
properties  on  both  sides  of  a  canopy  com¬ 
ponent  are  assumed  equal.  The  mul¬ 
tiplicative  separation  of  the  geometry  de¬ 
pendent  terms  from  the  energy  terms 
permits  the  possibility  of  precalculating 


the  geometric  matrix  S  for  a  wide  variety 
of  plant  canopies  as  a  means  of  char¬ 
acterizing  their  long-wave  thermal  behav¬ 
ior.  These  precalculated  matrices  may 
then  be  convolved  with  the  appropriate 
meteorological  driving  variables  as  re¬ 
quired  in  order  to  simulate  a  multitude  of 
target /background  scenarios.  In  addition, 
the  linear  form  of  Eq.  (2)  with  respect  to 
S  suggests  the  possibility  of  applying  lin¬ 
ear  filtering  theory  to  estimate  S  from 
measurements  of  canopy  temperature,  X. 

The  elements  of  the  Sj(  matrix  describe 
the  fraction  of  long-wave  emitted  flux 
from  a  source  layer,  i  (which  includes  the 
air  above  the  canopy  and  ground  layers  as 
well  as  the  three  vegetation  layers),  that  is 
intercepted  by  a  sink  layer  i.  This  flux 
must  escape  the  specific  source  layer  >' 
and  then  pass  unimpeded  through  all  in¬ 
termediate  layers  between  the  source  layer 
1  and  the  sink  layer  i  before  being  in¬ 
tercepted  by  the  canopy  elements  in 
layer  i. 

To  calculate  S(/,  we  integrate  over  all 
emitting  directions  0r,<j>r  the  total  flux 
that  escapes  a  source  layer  ;  and  is  in¬ 
tercepted  by  a  foliage  element  in  layer  i 
that  has  an  orientation  direction  de¬ 
scribed  by  foliage  inclination  angle  dk.  To 
calculate  the  total  flux  intercepted  in  layer 
i  from  a  source  layer  /,  we  then  sum  the 
energy  intercepted  by  each  foliage  in¬ 
clination  class  over  the  foliage  inclination 
angle  distribution  occurring  in  layer  i. 
Specifically, 

s„=  2  fikcltk,  do) 

k=l 

where 

flk  =  the  leaf  slope  distribution  for 
layer  i  =  1,2,3  and  foliage  in- 
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clinaticn  angle  dk  =  5°,  15°, 
...,85°. 

Ciik  =  fraction  of  emitted  flux  from  a 
source  layer  f  that  is  intercepted 
by  a  foliage  element  inclined  at 
angle  dk  within  layer  i. 

If  t  represents  a  unit  vector  in  the  direc¬ 
tion  of  an  emitting  source  element,  de¬ 
scribed  by  8r,<j>r,  and  <3,  a  unit  vector 
describing  the  orientation  of  an  absorbing 
sink  foliage  element,  then  the  amount  of 
flux  intercepted  by  this  foliage  element 
from  direction  t  is  proportional  to|d-f).  If 
we  let  CONTj(r  represent  the  fraction  of 
long-wave  flux  that  is  emitted  from  layer  j 
along  direction  f  and  is  finally  intercepted 
in  layer  t,  then,  the  total  flux  in¬ 
tercepted  by  the  foliage  element  emitted 
over  all  directions  f  in  layer  is  given  by 

c„,  =  /*/2/2i,|(9  .  fjcONT,  8r . 

Jo  Jo 

(11) 

Several  investigators  have  described  the 
calculation  procedure  for  estimating 
CONTj/r  for  various  theoretical  canopies, 
e.g.,  deWit  (1965)  and  Verhoef  and 
Bunnik  (1975).  Basically  the  calculations 
require  the  determination  of  the  probabil¬ 
ity  of  encountering  a  gap  (or  hit)  in 
traversing  a  vegetation  layer  which  is 


populated  by  foliage  elements  possessing 
a  leaf  slope  distribution  f(lt  and  a  leaf  or 
foliage  area  index.  We  have  generalized 
these  ideas  to  multiple  layers  (Oliver  and 
Smith,  1974).  The  specific  expressions  for 
these  weighting  coefficients  for  an  arbi¬ 
trary  source  direction  t  are  summarized  in 
Table  1.  In  this  table,  P0(i,  r)  represents 
the  probability  of  encountering  a  gap 
along  direction  t  in  traversing  layer  i. 
Note  that  the  probability  of  traversing 
half  a  layer  is  given  by  Pnl/2(i,  r)  and  that 
the  probability  of  encountering  emitting 
foliage  elements  in  a  layer  along  a  direc¬ 
tion  f  is  given  by  1  — P0(i,r).  For  exam¬ 
ple,  CONT13i.  represents  the  flux  that  is 
emitted  along  direction  f  from  source  layer 
3  (which  is  vegetation  canopy  layer  2) 
and  is  intercepted  by  sink  layer  1  (which 
is  vegetation  canopy  layer  1).  This  is  equal 
to  the  probability  of  having  emitting  ele¬ 
ments  in  canopy  layer  2:  i.e.,  the  proba¬ 
bility  of  a  hit  in  layer  2,  [1  — P0(2,  r)], 
times  the  probability  that  this  flux  en¬ 
counters  a  gap  in  traversing  to  the  mid¬ 
elements  of  absorbing  foliage  in  layer  1, 
P0l/2(1,  r).  It  is  thus  given  by 

CONT13r  =P01/2(l,r)[l  — P0(2,  r)] 

=  P01/2(l,r)-P01/2(l,r)P0(2,r) 

(12) 


TABLE  1  Expressions  lor  contribution  coefficients  CONTj(r  for  sink  layer  i,  so\irce  component  /,  and  arbitrary 
direction  6r.  P0(l,  r  )  =  probability  of  gap  in  layer  I  in  direction  $r. 


„  Sink  Layer 

Source  _ 

Layer  123 


1 

2 

3 

4 


P,  l/2(l,r) 

2[  1  —  Po /2(  l,  r )] 

P01/2{l,r)-Pn'/2(l.r)Pn<2,r) 
P<}/8(1. r  Wo(2.  r ) 

-  P<{ /2( L  r  Wri(2,  t )Pn(3. r ) 
Pn,/a(l.r)Pn(2.r)/'n(3.r) 


P(,<l.r)P,}/l(2.r> 

Po /a(2,  r)-P(}/a(2,  r)/o(l,  r) 

2(1  —  i’o /a(2,  r )] 
Pftl/i!(2,r)-P01/2(2.r)P11(3.r) 


P0(l,r)P0(2,r)P01/2(3.r) 
P«/2(3,  r  )P(,(2,  r ) 
-P01/2(3,r)P0(2.r)P„(l,r) 
P„1/2(3,r)-P01/2(3.r)P„(2,r) 
2[l-P„,/2(3.r)] 


5 


P„1/2(2,r)Pn(3,r) 


P(}/2(3,  r) 
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Model  Validation  Experiments 

Data  obtained  at  two  existing  research 
sites  were  utilized  to  evaluate  the  validity 
of  the  model.  One  site  is  located  in  a 
Douglas  fir  canopy  in  the  Cedar  River 
watershed  near  Seattle,  Washington.  The 
second  site,  at  the  Walker  Branch 
Watershed  near  Oak  Ridge,  Tennessee,  is 
typical  of  an  Appalachian  mixed  - 
deciduous  forest.  Both  research  sites  were 
being  used  for  ongoing  research  in  forest 
meteorology  and  had  extensive  instru¬ 
mentation  and  computerized  data 
acquisition  support. 

Cedar  River,  Washington  site 

The  Cedar  River,  Washington  study 
site  is  located  on  the  A.  E.  Thompson 
Research  Center  at  a  micrometeorological 
observatory  maintained  and  operated  by 
the  University  of  Washington  55  km 
southeast  of  Seattle,  Washington.  The 
average  elevation  is  approximately  215  m 
above  sea  level. 

The  dominant,  naturally  regenerated 
stand  of  Douglas  fir  [  Pseudotsuga 
menziesii  (Mirb.)  Franco]  was  approxi¬ 
mately  41  years  old  with  an  average  tree 
spacing  of  5.8  m.  Average  height  of  the 
Douglas  fir  stand  was  about  28  m  with  an 
average  leaf-area  index  (LAI)  of  ap¬ 
proximately  7.8.  Ground  cover  consisted 
of  fern,  salal,  huckleberry,  mosses,  and 
litter.  Soil  at  the  site  consisted  of  Bames- 
ton  gravelly  loamy  sand  originating  from 
glacial  outwash  (Jensen,  1976), 

Located  at  this  site  was  a  28-m  tall 
Douglas  fir  tree  contained  in  a  lysimeter 
(Fritschen  et  al.,  1973).  The  site  adjacent 
to  this  tree  was  instrumented  to  provide 
data  for  evapotranspiration  studies.  These 
data  included  wet  and  dry  bulb  tempera¬ 
ture  profiles,  soil  temperatures,  global 
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short-wave  radiation,  precipitation,  and 
wind  speed  and  direction.  In  addition, 
needle  surface  temperatures  were  moni¬ 
tored  at  several  points  around  the  lysime¬ 
ter  tree  near  the  top  and  center  of  the 
canopy. 

Walker  Branch,  Tennessee  site 

The  Walker  Branch  study  site  is  located 
near  the  Walker  Branch  Watershed  re¬ 
search  facility  on  the  U.S.  Department  of 
Energy  Reservation  near  Oak  Ridge,  Ten¬ 
nessee.  This  research  area  is  situated  on  a 
ridge  top  about  70  m  above  the  valley 
floor  at  an  elevation  of  335  m  above  mean 
sea  level. 

The  area  is  representative  of  an  Appa¬ 
lachian  deciduous  forest  (Hutchison, 
1977).  The  species  composition  of  the 
stand  is  dominated  by  various  species  of 
oak  and  hickory.  The  average  height  of 
the  codominant  trees  is  about  21.5  m 
with  lower  limit  of  the  live  crown  being 
15  m  above  the  ground.  Basal  area  was 
approximately  26  m2  ha-1.  Understory 
growth  is  abundant  and  the  groimd  is 
covered  by  a  shallow  accumulation  of 
litter.  Hutchison  (1977)  gives  a  detailed 
description  of  the  site  and  data  available 
at  the  research  facility. 

Modeling  Input  Data 

The  data  collected  at  the  two  sites  in¬ 
clude  foliage  and  background  optical 
parameters,  geometry  characterization 
measurements,  and  environmental  mea¬ 
surements.  This  section  describes  the  data 
required  for  our  models  and  the  tech¬ 
niques  or  sources  used  to  acquire  it. 

Foliage  geometry 

The  procedure  for  determining  foliage 
inclination  angles  for  the  Douglas  fir 
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canopy  included  acquiring  high  contrast 
black  and  white  slide  photography  of 
canopy  silhouettes.  For  the  purposes  of 
our  modeling,  the  canopies  were  parti¬ 
tioned  into  three  layers  of  equal  height. 
High  contrast  slides  were  used  as  input  to 
a  laser  diffractometer  and  the  diffraction 
patterns  then  optically  sampled  (Kimes  et 
al.,  1979).  Separate  branch  and  foliage 
measurements  were  combined  to  provide 
the  inclination  angle  distributions  for  each 
layer. 

Leaf  inclination  distributions  for  the 
oak-hickory  canopy  were  sampled  in  situ 
for  several  leaves  at  1-m  intervals 
throughout  the  stand.  The  recorded  dis¬ 
tributions  were  summed  and  averaged 
over  the  appropriate  layer-height  intervals 
to  provide  the  three-layer  leaf-inclination 
angle  distributions.  The  leaf-angle  distri- 
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buttons  for  the  two  canopies  are  com¬ 
pared  on  a  layer  by  layer  basis  in  Fig.  1. 

Leaf-area  index  (LAI)  is  defined  as  the 
total  one-sided  leaf  area  occupying  the 
horizontally  projected  area  of  the  canopy. 
LAIs  for  the  Douglas  fir  canopy  were 
derived  from  measurements  reported  by 
Kinerson  and  Fritschen  (1971).  In  this 
paper,  graphs  of  canopy  height  z(m) 
versus  surface  area  density  F{z),(  rrPnr*) 
for  nine  sample  plots  are  given.  Integrat¬ 
ing  F(z)  over  height  gives  the  needle- 
surface  area  index  (NSAI)  for  a  particular 
height  increment  ( dz ).  For  our  modeling 
purposes,  LAI  values  were  determined  by 
dividing  NSAI  for  each  layer  by  two. 
Data  collected  at  the  site  since  these  mea¬ 
surements  were  made  indicate  no  sub¬ 
stantial  change  of  LAI  since  1971.  Values 
of  LAI  for  the  oak-hickory  canopy  layers 


FIGURE  I.  Comparative  plot*  of  foliage  inclination  angle  vs.  cumulative  frequency  for  the  three  layer  Douglas  fir 
and  oak-hickory  canopies.  A  =  Layer  1.  B  =  Layer  2,  C  =  Layer  3. 


[•Te] 


LEAF  INCLINATION  ANGLE  (degrees) 
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were  derived  from  a  graph  of  cumulative 
LAI  versus  height  through  the  canopy. 
The  graphs  were  generated  by  ATDL 
personnel  from  direct  measurements.  Ta¬ 
ble  2  lists  the  LAIs  and  midlayer  heights 
used  to  model  the  Douglas  fir  and  oak- 
hickory  canopies  and  the  S  matrices  for 
each  canopy. 

Short-wave  absorption  coefficients 

The  absorption  of  global  short-wave 
radiation  by  canopy  layers  is  an  im¬ 
portant  component  in  the  daytime  energy 
budget.  The  radiation  absorbed  is  a  func¬ 
tion  of  the  global  short-wave  energy  avail¬ 
able,  which  is  a  measured  input  parame¬ 
ter  to  the  model;  and  the  short-wave  in¬ 
terception  coefficients  for  the  canopy  sys¬ 
tem,  which  in  our  case  were  estimated  by 
a  separate  multiple  scattering  absorptive 
Monte  Carlo  model  (Kimes  and  Smith, 
1980).  Strictly  speaking  these  short-wave 
absorption  coefficients  for  the  canopy  and 
ground  layers  are  a  function  of  sun  angle, 
a  result  borne  out  by  the  Monte  Carlo 
analyses.  However,  for  the  thermal  model 
reported  here  our  objective  was  to  sim¬ 
plify  the  required  analyses  in  order  to 


facilitate  the  analysis  of  a  multitude  of 
canopy  situations  and  reduce  required  in¬ 
put  data.  Extensive  analyses  with  the 
complete  treatment  of  short-wave  energy 
absorption  would  be  expensive  and  previ¬ 
ous  analyses  indicate  that  the  absorption 
coefficients  were  relatively  stable  within 
15-20%  for  sun  angles  ranging  from  nadir 
to  zenith  angles  of  45°.  It  was  felt  that  an 
average  absorption  value  calculated  for 
each  canopy  for  these  sun  angles  gener¬ 
ally  would  be  reasonable.  At  large  zenith 
solar  angles  the  short-wave  absorption 
coefficient  becomes  highly  nonlinear.  In 
our  treatment  of  this  parameter  we  tend 
to  underestimate  this  component  of  the 
energy  source  term  significantly  at  early 
morning  and  late  evening,  but  the  magni¬ 
tude  of  the  insolation  is  relatively  small  at 
these  hours. 

In  order  to  estimate  the  short-wave  flux 
absorption  it  is  necessary  to  have  esti¬ 
mates  of  the  canopy  and  ground  optical 
scattering  properties,  i.e.,  reflectance  and 
transmittance  of  foliage  elements  and 
ground  layer. 

Canopy  element  transmittance  values 
were  directly  measured  at  both  sites  as 
well  as  average  background  reflectance 


TABLE  2  Canopy  layer  heights,  LAIs  anil  S  matrices  for  the  Douglas  fir  and 
oak-hickory  canopies  modeled  in  this  study. 


Mid-layer  height  (m) 

Leaf-area  index 

Layer  1 

Layer  2 

Layer  3 

Layer  1 

Layer  2 

Layer  3 

To  (i)  Layer 

Douglas  fir 

23,3 

14,0 

4.7 

1.5 

5.3 

1.0 

Oak 

-HICKORY 

18.3 

11.0 

3.7 

3.4 

0.8 

0.4 

To  (i)  Layer 

S  Matrix 

l 

2 

3 

1 

2 

3 

From  (f) 

Sky 

0,2722 

0.0006 

0.0000 

0.1595 

0,0281 

0.0201 

Layer  1 

1.4484 

0.0048 

0.0000 

1.6741 

0.7914 

0.5441 

Layer  2 

0.2722 

1.9820 

0.2946 

0,0470 

0.3539 

0.2574 

Layer  3 

0.0000 

0.0047 

1.4035 

0.0338 

0.2589 

0,3496 

Ground 

0.0000 

0,0007 

0.2946 

0.0788 

0.5607 

0.8217 
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values  as  described  In  the  report  (Smith 
et  ah,  1981).  However,  it  was  not  practi¬ 
cal  to  obtain  measurements  of  the  canopy 
element  reflectances,  particularly  for  the 
Douglas  fir  needles.  Rather,  literature  re¬ 
flectance  values  for  both  old  and  new 
Douglas  fir  were  obtained  from  Jarvis 
et  al.  (1976).  Similarly,  measurements  by 
Colwell  (1969)  were  averaged  for  the 
oak -hickory  leaves.  In  both  cases  the  ac¬ 
tual  site  measured  transmittance  values 
were  used  to  ensure  that  at  least  physi¬ 
cally  reasonable  reflectance  estimates 
were  obtained. 

Stomatal  resistance 

The  resistance  of  the  leaf  to  water  vapor 
diffusion  depends  on  many  environmental 
factors.  Lead  stomates  open  and  close  in 
response  to  microclimatic  and  soil  condi¬ 
tions  and  regulate  the  cooling  of  the  plant 
through  evapotranspiration.  This  parame¬ 
ter  is  difficult  to  measure  and  highly  vari¬ 
able.  For  our  modeling  purposes  average 
values  were  used  as  constants.  The  value 
for  Douglas  fir  was  set  at  0.10  min/cm. 
Stomatal  resistances  measured  for  the 
oak-hickory  canopy  ranged  from  0.04  to 
0.07  min/cm  for  sun  leaves.  The  upper 
value  was  selected  for  use  in  the  decidu¬ 
ous  canopy  simulations.  Stomatal  resis¬ 
tance  was  set  to  infinity  during  nighttime 
hours  for  both  canopies.  Model  estimates 
of  leaf  temperatures  are  not  very  sensitive 
to  stomatal  resistance  at  these  values  and 
under  the  moderate  environmental  condi¬ 
tions  encountered  during  the  data  collec¬ 
tion  (Smith  et  al.,  1981). 

Emissivity  and  absorptivity 

The  ability  of  a  canopy  element  to  emit 
and  absorb  long-wave  radiation  is  ex¬ 
pressed  by  the  emissivity  and  absorptivity 
coefficients  specified  in  the  model.  Emis¬ 


sivity  (e,)  and  absorptivity  (a,)  were  set 
to  1.0  for  all  three  canopy  layers.  Emissiv¬ 
ity  of  the  ground  (eK)  was  also  set  to  1.0. 
Emissivity  of  the  air  (ett)  was  calculated 
as  a  function  of  air  temperature  by  an 
empirical  relationship  (Hudson,  1969). 

Canopy  temperature  measurements 

Since  the  purpose  of  the  experiments 
was  to  collect  data  sets  for  validation  of 
the  thermal  ukkN,  actual  canopy  foliage 
temperature  measurements  were  re¬ 
quired. 

The  experimental  setup  at  the  Cedar 
River  site  included  temperature  measure¬ 
ments  for  a  number  of  individual  Douglas 
fir  needles.  The  temperature  sensors  were 
located  around  the  lysimeter  tree  at 
heights  from  20  to  26  m.  The  measure¬ 
ments  at  a  given  height  were  averaged  to 
give  an  average  layer  temperature.  The 
26-m  measurement  was  assumed  to  repre¬ 
sent  the  average  canopy  temperature  for 
the  top  layer  (Layer  1).  The  20-m  mea¬ 
surement  was  assumed  to  approximate 
the  middle  layer  (Layer  2)  although  its 
location  is  closer  to  the  boundary  be¬ 
tween  Layer  1  and  Layer  2. 

No  individual  leaf  temperature  mea¬ 
surements  were  available  at  the  Walker 
Branch  site,  so  a  portable  thermal  radiom¬ 
eter1  was  used  to  monitor  the  canopy 
temperature  throughout  a  24-hr  period. 
The  procedure  was  to  position  the  instru¬ 
ment  upward  from  the  ground  at  the 
canopy  and  slowly  move  it  until  the  maxi¬ 
mum  temperature  was  recorded.  This  was 
done  to  minimize  errors  due  to  the  pres¬ 
ence  of  sky  or  clouds  in  the  field  of  view. 
However,  the  inferred  measured  tempera¬ 
ture  represents  a  value  integrated  over 
the  entire  depth  of  the  canopy  and 

5  Barnes  Insta-Therm,  Barnes  Engineering  Corporation 
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weighted  most  heavily  toward  the  bottom  peratures.  All  measurements  were  either 
of  the  canopy.  It  is,  thus,  not  a  precise  instantaneous  or  short  time  interval  aver- 
measurement.  ages. 


Environmental  input  parameters 

In  addition  to  the  geometrical,  optical 
and  thermal  parameters  discussed  above, 
a  set  of  dynamic  variables  characterizing 
the  microclimate  of  the  target  are  re¬ 
quired  to  drive  the  thermal  model.  These 
parameters  consist  of  air  temperature 
above  the  canopy,  ground  surface  tem¬ 
perature,  wind  speed  at  the  top  of  the 
canopy,  relative  humidity,  and  global 
short-wave  radiation. 

Environmental  data  were  provided 
from  the  automated  recording  systems  at 
the  two  sites.  Air  and  ground  tempera¬ 
tures  and  global  short-wave  radiation  were 
measured  directly.  Relative  humidity  was 
determined  from  wet  and  dry  bulb  tem- 


Results 

The  data  collected  for  the  coniferous 
and  deciduous  canopies  provided  a  good 
means  of  testing  the  thermal  model  under 
these  diverse  conditions.  Three  layer 
canopy  temperature  simulations  were 
made  over  a  48-hr  period  with  both  data 
sets  and  the  results  were  compared  with 
measured  temperatures. 

Douglas  fir  canopy 

The  thermal  model  was  run  with  en¬ 
vironmental  data  acquired  over  the  48-hr 
period  of  4-5  August  1979. 

The  Layer  1  simulated  temperatures 
followed  the  trend  of  air  temperature 


4  AUGUST  1979  I  5  AUGUST  1979 

FIGURE  2.  Layer  1  predicted  temperatures  plotted  with  average  temperatures 
measured  at  the  26-m  level  in  the  Douglas  fir  canopy. 
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TIME  (hours) 

4  AUGUST  1979  I  5  AUGUST  1979 

FIGURE  3.  La/er  2  £r«!(ctc'i  temperature  plotted  with  average  temperatures 
measured  at  the  20- tn  level  in  the  Douglas  lir  canopy. 


throughout  the  48-hr  period.  Comparisons 
of  measured  and  predicted  needle  tem¬ 
peratures  are  presented  as  Fig.  2  and  3 
for  Layers  1  and  2,  respectively. 

The  Layer  1  predicted  temperatures 
varied  from  measured  by  a  maximum  of 
3°C .  These  deviations  were  observed  dur¬ 
ing  the  daylight  hours  under  hazy  skies. 
Nighttime  predictions  deviated  from 
measured  by  2°C  or  less  with  the  maxi¬ 
mum  deviations  occurring  under  condi¬ 
tions  of  fog.  This  leads  us  to  conclude 
that  the  thermal  model  may  be  most  valid 
for  days  with  primarily  direct  solar  radia¬ 
tion  and  clear  nights  where  radiative  cool¬ 
ing  is  occurring. 


validate  the  thermal  model  for  a  decidu¬ 
ous  oak-hickory  canopy.  Nighttime  simu¬ 
lations  were  nearly  equal  to  air  tempera¬ 
ture  while  daytime  predictions  varied  by 
a  maximum  of  2°C  over,  air  temperature. 

Measured  temperatures  were  compared 
to  predicted  results  for  Layer  2  and  are 
shown  in  Fig.  4.  The  agreement  between 
model  and  measured  temperatures  is  quite 
good.  However,  as  discussed  earlier  an 
indirect  measure  of  canopy  effective 
radiant  temperatures  was  made.  The 
largest  deviation  (3°C)  occurs  in  the 
afternoon  whereas  morning  and  night¬ 
time  predictions  vary  only  1°C  or  less. 

Summary 


Oak-hickory  canopy 


A  simplified  thermal  canopy  model 


Environmental  data  acquired  at  the  which  treats  the  radiative  flux  transfers  in 
Walker  Branch  site  for  the  48-hr  period  some  detail  and  permits  the  incorporation 
from  18-19  August  1979  were  used  to  of  alternative  expressions  for  other  energy 
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TIME  (hours) 

IS  AUGUST  1979  I  19  AUGUST  1979 


FIGURE  4.  Layer  2  predicted  temperature  plotted  measured  average  tempera¬ 
ture  of  oak-hickory  canopy.  Measurements  were  made  from  the  ground  with  a 
thermal  IR  radiometer. 

components  has  been  described  and  measured  observations.  This  latter  sugges- 
shown  to  give  reasonable  results  for  two  tion  is  currently  being  intensely  investi- 
forest  canopy  situations.  The  model  gated  by  the  authors, 
utilized  detailed  canopy  structure  char-  The  results  of  the  model-experiment 
acteristics  but  in  contrast  to  an  earlier  but  comparisons  indicate  that  the  model  pro- 
more  comprehensive  form  of  the  model  vides  reasonable  estimates  of  actual  tem- 
(Kimes  et  al.,  1981)  factors  the  long-wave  peratures  for  nighttime  periods  tc  within 
energy  source  terms  into  a  product  of  two  2°C  for  both  canopies  studied.  Daytime 
terms,  one  dependent  on  thermal  proper-  simulations  generally  deviated  from  mea- 
ties  only  of  the  canopy  (the  Boltzmann  sured  temperatures  because,  perhaps,  of 
source  term)  and  the  other  on  canopy  the  simplifications  assumed  for  the  short- 
geometric  characteristics  only.  This  fac-  wave  flux  absorption  coefficients  and  the 
torization  has  two  significant  applica-  treatment  of  stomatal  resistance.  The  re- 
tions.  First,  as  was  done  here,  geometrical  suits  indicate  that  the  model  may  not 
long-wave  flux  transfer  matrices  may  be  adequately  account  for  some  of  the  en- 
precalculated  for  a  variety  of  targets  and  ergy  transfers  under  more  extreme  en- 
subsequently  convolved  with  local  vironmental  conditions, 
meteorological  scenarios.  Secondly,  the  Finally,  it  should  be  noted  that  both  in 
linear  form  of  the  factorization  suggests  measurement  situations,  the  model  pre- 
the  use  of  linear  filtering  theory  to  esti-  dictions,  air  temperature,  and  canopy 
mate  these  flux  transfer  matrices  given  temperature  measurements  all  agree  quite 
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closely.  The  thermal  model,  however,  pro¬ 
vides  a  convenient  organization  of  the 
energy  flow  processes  in  a  self-consistent 
manner  which  relates  measured  or  pre¬ 
dicted  canopy  temperatures  to  intrinsic 
canopy  parameters.  It  thus  permits  the 
potential  inference  of  canopy  characteris¬ 
tics  from  measurements  as  suggested 
above. 

The  research  described  in  this  paper 
was  supported  by  the  U.  S.  Army  Corps  of 
Engineers  Waterways  Experiment  Station 
under  Contract  DACW  39-77-C-0073  and, 
in  part,  by  the  C/.S.  Army  Research 
Office  —  Durham  under  Grant  No. 
DAAG29 — 79-C-0199.  We  wish  to  ac¬ 
knowledge  the  helpful  assistance  of  Al¬ 
fonso  Vasques  of  the  WntewMys  Experi¬ 
ment  Station  with  data  collection  experi¬ 
ments. 

References 

Alderfer,  R.  G.,  and  Gates,  D.  M.  (1971), 
Energy  exchange  in  plant  canopies,  Ecol¬ 
ogy  52(5):855- 861. 

Burden,  R.  L.,  Faires,  J.  D.,  ar.d  Reynolds, 
A.  C.  (1978),  Numerical  Analysis,  Prindle, 
Weber  and  Schmidt,  Boston,  MA. 

Colwell,  E.  E.  (1969),  Seasonal  change  in 
foliar  reflectance  of  five  broadleaved  forest 
tree  species,  Ph.D.  Thesis,  University  of 
Michigan,  Ann  Arbor,  112  p. 

Deardorff,  J.  W.  (1978),  Efficient  prediction 
of  ground  surface  temperature  with  inclu¬ 
sion  of  a  layer  of  vegetation,  J.  Geophys. 
Res.  (83):  1889- 1904. 

deWit,  C.  T.  (1965),  Photosynthesis  of  leaf 
canopies.  Agr.  Res.  Pap.  663,  Wageningen, 
Netherlands,  pp.  1-57. 

Fritschen,  L.  J.,  Cox,  L.,  and  Kinerson,  R. 
(1973),  A  28-meter  Douglas-fir  in  a  weigh¬ 
ing  lysimeter,  Forest  Sci.  19:256-261. 


Gates,  D.  M.  (1968),  Energy  Exchange  in  the 
Biosphere,  Harper  and  Row,  Inc.,  New 
York,  131  p. 

Goudriaan,  J.  (1977),  Crop  micrometeorology: 
a  simulation  study,  Simulation  Monographs 
PUDOC,  Centre  for  Agricultural  Publish¬ 
ing  and  Documentation,  Wageningen,  The 
Netherlands,  249  p. 

Heilman,  J.  L.,  Kanernasu,  E.  T.,  and  Rosen- 
berg,  J.  J.  (1976),  Thermal  scanner  mea¬ 
surement  of  canopy  temperatures  to  esti¬ 
mate  evapotranspiration,  Remote  Sens.  En¬ 
viron.  5(2):  137- 145. 

Hudson,  R.  D.,  Jr.  (1969),  Infrared  System 
Engineering,  Wiley,  New  York. 

Hutchison,  B.  A.  (1977),  Atmospheric  turbu¬ 
lence  and  diffusion  laboratory  deciduous 
forest  meteorology  research  program,  an 
overview,  ATDL,  NOAA,  Cont.  File  77/1, 
lip. 

Jarvis,  P.  G.,  James,  G.  B.,  and  Landsberg,  J. 
J.  (1976),  in:  Vegetation  and  the  Atmo¬ 
sphere,  Vol.  2  (J.  L.  Monteith,  Ed.), 
Academic  Press,  New  York,  439  p. 

Jensen,  E.  C.  (1976),  The  crown  structure  of  a 
single  codominant  Douglas-fir,  M.S.  Thesis, 
University  of  Washington,  Seattle, 
Washington,  83  p. 

Kimes,  D.  S.,  Ranson,  K.  J.,  Kirchner,  J.  A., 
and  Smith,  J.  A.  (1978),  Modeling  descrip¬ 
tors  and  terrain  modules.  Final  Report  un¬ 
der  contract  DACW  29-77-C-0073,  En¬ 
vironmental  Laboratory,  U.S.  Army  En¬ 
gineers  Waterways  Experiment  Station, 
Vicksburg,  Mississippi  125  p. 

Kimes,  D.  S.,  Smith,  J.  A.,  and  Ranson,  K.  J. 
(1979),  Terrain  feature  canopy  modeling. 
Final  Report  under  U.S.  Army  Research 
Office  Grant  DAAG29-78-0045,  Colorado 
State  University,  Fort  Collins,  CO. 

Kimes,  D.  S.,  Smith,  J.  A.  and  Berry,  K.  J. 
(1979),  Extension  of  the  optical  diffraction 
analysis  technique  for  estimating  forest 
canopy  geometry,  Aust.  ].  Bot.  27:575-588. 

Kimes,  D.  S.,  and  Smith,  J.  A.  (1980),  Sinru- 


326 


J.  A.  SMITH  ET  AL. 


lation  of  solar  radiation  absorption  in  vege¬ 
tation  canopies,  Appl.  Opt.  19(16):2801- 
2811. 

Kimes,  D.  S.,  Ranson,  K.  J.,  and  Smith,  J.  A. 
(1980),  A  Monte-Carlo  calculation  of  the 
effects  of  canopy  geometry  on  PhAR  ab¬ 
sorption,  Photosynthetica  14(1):55— 64. 

Kimes,  D.  S.,  Smith,  J.  A.  and  Link,  L-  E. 
(1981),  A  thermal  IR  exitance  model  of  a 
plant  canopy.  Appl.  Opt.  20(4);623-632. 

Kinerson,  R.  Jr.,  and  Fritschen,  L.  J.  (1971), 
Modeling  a  coniferous  forest  canopy,  Agr. 
Meteor.  8:439-445, 

Norman,  J.  (1979),  Modeling  the  cor  plete 
crop  canopy.  In:  Modification  of  the  En¬ 
ergy  Environment  of  Plant  (B.  J,  Barfield 
and  J  F.  Gerber,  Eds.),  ASAE  Monograph 
Number  2,  Amer.  Soc.  of  Ag.  Engs.,  St. 
Joseph,  MI,  539  p. 

Oliver,  R.  E.,  and  Smith,  J.  A.  (1974),  A 
stochastic  canopy  model  of  diurnal  reflec¬ 
tance,  Final  Report,  U.S.  Army  Research 
Office,  Durham,  NC. 

Reginato,  R.  J.,  Idso,  S.  B.,  Vedder,  J.  F., 
Jackson,  R.  D.,  Blanchard,  M.  B.,  and 
Goettelman,  R.  (1970),  Soil  water  content 
and  evaporation  obtained  from  ground- 
based  and  remote  measurements,  J.  Geo- 
phys.  Res.  81(9):  1617- 1620. 

Soer,  F.  J.  R.  (1980),  Estimation  of  regional 
evapotranspiration  and  soil  moisture  condi¬ 


tions  using  remotely  sensed  crop  surface 
temperatures.  Remote  Sens.  Environ.  9:27- 
45. 

Sorenson,  H.  W.  (1966),  Kalman  Filtering 
Techniques,  Advances  in  Control  Systems, 
Vol.  3,  219  (C.  T.  Leondes,  Ed.),  Academic 
Press,  New  York. 

Smith,  J.  A.,  Ranson,  K.  J.,  Nguyen,  D.,  and 
Link,  L.  E.  (1981),  Thermal  vegetation 
canopy  model  studies,  Final  Report  under 
contract  DACW39-77-C-0073,  Environ¬ 
mental  Laboratory,  U.S.  Army  Engineers, 
Waterways  Experiment  Station,  Vicksburg, 
MS,  297  p. 

Tibbals,  E.C.,  Carr,  E.K.,  Gates,  D.M.,  Kreith, 
F.  (1964),  Radiation  and  convection  in  con¬ 
ifers,  Amer.  J.  Bot.  51(5):529-538. 

Verhoef,  W.,  and  Bunnik,  N.  J.  J.  (1975),  A 
model  study  on  the  relations  between  crop 
characteristics  and  canopy  spectral  reflec¬ 
tance,  NIWARS  publication  No.  33,  3 
Knaalweg  Delft,  the  Netherlands,  89  p. 

Wiebelt,  J.  A.,  and  Henderson,  J.  B.  (1977), 
Techniques  and  analysis  of  thermal  in¬ 
frared  camouflage  in  foliated  backgrounds, 
Final  Report  under  contract  DAAG-76-C- 
0134,  U.S.  Army  Mobility  Equipment  Re¬ 
search  and  Development  Command,  Fort 
Belvoir,  VA,  63  p. 

Received  14  April  1980;  revised  16  January  1981. 


